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We review five mechanisms for forming brown dwarfs: (i) turbulent fragmentation of 
molecular clouds, producing very low-mass prestellar cores by shock compression; (ii) 
collapse and fragmentation of more massive prestellar cores; (iii) disc fragmentation; (iv) 
premature ejection of protostellar embryos from their natal cores; and (v) photo-erosion of 
pre-existing cores overrun by HII regions. These mechanisms are not mutually exclusive. 
Their relative importance probably depends on environment, and should be judged by their 
ability to reproduce the brown-dwarf IMF, the distribution and kinematics of newly formed 
brown dwarfs, the binary statistics of brown dwarfs, the ability of brown dwarfs to retain discs, 
and hence their ability to sustain accretion and outflows. This will require more sophisticated 
numerical modelling than is presently possible, in particular more realistic initial conditions 
and more realistic treatments of radiation transport, angular momentum transport and magnetic 
fields. We discuss the minimum mass for brown dwarfs, and how brown dwarfs should be 
distinguished from planets. 



1. INTRODUCTION 

The existence of brown dwarfs was first proposed on the- 
oretical grounds by Kumar (1963) and Hayashi and Nakano 
(1963). However, more than three decades then passed be- 
fore brown dwarfs were observed unambiguously (Rebolo 
et al., 1995; Nakajima et al., 1995; Oppenheimer et al., 
1995). Brown dwarfs are now observed routinely, and it is 
therefore appropriate to ask how brown dwarfs form, and 
in particular to ascertain (a) whether brown dwarfs form 
in the same way as H-burning stars, and (b) whether there 
is a clear distinction between the mechanisms that produce 
brown dwarfs and those that produce planets. 

In §2 we argue that the mechanisms forming brown 
dwarfs are no different from those forming low-mass H- 
burning stars, on the grounds that the statistical proper- 
ties of brown dwarfs (mass function, clustering properties, 
kinematics, binary statistics, accretion rates, etc.) appear 
to form a smooth continuum with those of low-mass H- 
burning stars. Understanding how brown dwarfs form is 
therefore the key to understanding what determines the min- 
imum mass for star formation. In §3 we review the basic 



physics of star formation, as it applies to brown dwarfs, 
and derive key analytic results, in particular the minimum 
mass for opacity-limited fragmentation in various different 
formation scenarios. We have deliberately assembled most 
of the mathematical analysis in this one section. In §4 to 
§8 we consider five different mechanisms that may be in- 
volved in the formation of brown dwarfs. §4 explores the 
possibility that turbulent fragmentation of molecular clouds 
produces prestellar cores of such low mass that inevitably 
they collapse to form brown dwarfs. §5 considers more 
massive prestellar cores, and the possibility that they frag- 
ment dynamically as they collapse, thereby spawning pro- 
tostars with a range of masses. The collapse of such cores 
ceases when the gas starts to heat up (due to adiabatic com- 
pression) and/or when rotation becomes important. Un- 
fortunately, the interplay of dynamics and adiabatic com- 
pression, whilst likely to play a critical role (e.g. Boss 
et al, 2000), is hard to analyse. On the other hand ro- 
tational effects can be analysed systematically, and there- 
fore §6 explores the formation of brown dwarfs by gravita- 
tional instabilities in discs, considering first isolated relaxed 
discs, then unrelaxed discs, and finally interacting discs in 
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clusters. §7 reviews the process of competitive accretion, 
which determines how an ensemble of protostellar embryos 
evolves to populate higher masses; and then considers the 
iV-body processes which (a) may eject brown dwarfs and 
low-mass stars from their natal cores, thereby terminating 
accretion and effectively capping their masses, and (b) may 
influence the binary statistics and clustering properties of 
brown dwarfs. §8 considers the formation of brown dwarfs 
by photo-erosion of pre-existing cores which are overrun 
by HII regions. §9 presents numerical simulations of the 
birth of a whole star cluster in a large protocluster core; in 
these simulations, the mechanisms of §5 (collapse and frag- 
mentation), §6 (disc fragmentation) and §7 (premature ejec- 
tion) all occur simultaneously and in tandem, and the col- 
lective properties of brown dwarfs formed can be extracted 
for comparison with observation. In §10 we summarise the 
principal conclusions of this review. 

Comprehensive discussions of the observational proper- 
ties of brown dwarfs, and how they may constrain forma- 
tion mechanisms, are contained in the chapters by Luhman 
et al. (which deals with the entire observational picture) and 
Burgasser et al. (which deals specifically with the issue of 
multiplicity). Therefore we offer only a brief summary of 
the observations in §2. Likewise, a critique of simulations 
of core fragmentation is given in the chapter by Goodwin 
et ah; simulations of disc fragmentation are discussed and 
compared in detail in the chapter by Durisen et al. ; and the 
origin of the IMF is treated in the chapter by Bonnell et al, 
so we have limited our consideration to those aspects which 
pertain specifically to the origin of brown dwarfs. 

2. EVIDENCE THAT BROWN DWARFS FORM 
LIKE LOW-MASS H-BURNING STARS 

We shall assume that brown dwarfs form in the same 
way as H-burning stars, i.e. on a dynamical timescale, by 
gravitational instability, and with initially uniform elemen- 
tal composition (reflecting the composition of the interstel- 
lar matter out of which they form). (By implication, we 
distinguish brown dwarfs from planets, a term which we 
will here reserve for objects which form on a much longer 
timescale, by the amalgamation of a rocky core and - if cir- 
cumstances allow - the subsequent accretion of a gaseous 
envelope. This results in planets having an initially fraction- 
ated elemental composition with an overall deficit of light 
elements.) If this is the correct way to view the formation 
of brown dwarfs, then brown dwarfs should not be distin- 
guished from stars. Many stars fail to burn helium, and most 
fail to burn carbon, without forfeiting the right to be called 
stars. The reason for categorising brown dwarfs as stars 
is that the statistical properties of brown dwarfs appear to 
form a continuum with those of low-mass H-burning stars 
(and not with those of high-mass planets). 

2.1. The IMF, clustering statistics, velocity dispersion 

The initial mass function (IMF) is apparently continu- 
ous across the hydrogen-burning limit at ~ 0.075 M , in 



the Trapezium cluster (Slesnick et al., 2004), in the a On- 
onis cluster (Bejar et al., 2001), in Taurus (Luhman, 2004), 
in IC348 (Luhman et al, 2003), in the Pleiades (Moraux 
et al, 2003), in the field stars of the disc (Chabrier, 2003) 
and even possibly in the halo (Burgasser et al, 2003b). If 
the IMF is fitted by a power law across the H-burning limit, 
dJV/dM oc M~ a , estimates of a fall in the range 0.4 to 0.8 
(e.g. Moraux et al, 2003). The IMF appears to extend down 
to a few Jupiter masses (e.g. Zapatero Osorio et al, 2002; 
McCaughrean et al, 2002; Lucas et al, 2005). The conti- 
nuity of the IMF across the H-burning limit is not surpris- 
ing, since the processes which determine the mass of a low- 
mass star are presumed to occur at relatively low densities 
(p < 10~ 8 g cm~ 3 ) and temperatures (T < 2000 K), long 
before the material involved knows whether it will reach 
sufficiently high density (p > 1 g cm~ 3 ) to be supported in 
perpetuity by electron degeneracy pressure before or after it 
reaches sufficiently high temperature (T > 10 7 K) to burn 
hydrogen. 

In the Trapezium cluster, and in Taurus, brown dwarfs 
appear to be homogeneously mixed with H-burning stars 
(Lucas and Roche, 2000; Briceno et al, 2002; Luhman, 
2004), and in Chamaeleon I their kinematics are also essen- 
tially indistinguishable (Joergens, 2006b). Although they 
have been searched for - as possible signatures of forma- 
tion by ejection - neither a greater velocity dispersion of 
brown dwarfs in very young clusters, nor a diaspora of 
brown dwarfs around older clusters, has been found. 

2.2. Binary statistics 

Multiple systems involving brown dwarf secondaries can 
be categorised based on whether the primary is a Sun-like 
star or another brown dwarf. 

Sun-like primaries. Sun-like stars seldom have brown- 
dwarf companions. At close separations (< 5AU), 
the frequency of companions with masses in the range 
0.01 to 0.1 M Q is - 0.5% (Marcy and Butler, 2000); 
this figure is known rather accurately due to the numer- 
ous Doppler surveys aimed at detecting extrasolar planets. 
Since the frequency of companions outside this mass range 
(both exoplanets below, and H-burning stars above) is much 
higher, the paucity of brown-dwarf companions is termed 
the 'brown dwarf desert' . However, the chapter by Luh- 
man et al. points out that the paucity of close brown-dwarf 
companions to Sun-like stars may simply reflect the overall 
paucity of brown dwarfs. At larger separations (> 100 AU), 
only about 20 brown-dwarf companions to Sun-like stars 
have been found to date, indicating that these systems are 
also rare (frequency ~ 1 to 2 %), though this estimate is 
based on more limited statistics (e.g. Gizis et al, 2001; 
McCarthy and Zuckerman, 2004). 

Brown-dwarf primaries. BD-BD binary systems (and 
binary systems with very low-mass H-burning primaries) 
have an observed multiplicity of ~ 10 — 20% for separa- 
tions greater than ~ 2 AU (Bouy et al, 2003; Burgasser et 
al, 2003a; Close et al, 2003; Gizis et al, 2003), and all but 
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five of the ~ 75 known BD-BD binaries have separations 
less than ~ 20 AU (see Table 1 in the chapter by Burgasser 
et al). Below ~ 2 AU, the multiplicity is unknown because 
most surveys to date have been imaging surveys that cannot 
resolve close systems. Only a few spectroscopic BD-BD bi- 
naries have been discovered (Basri and Martin, 1999; Stas- 
sun et al., 2006; Kenyon et al, 2005). Some authors spec- 
ulate that the overall multiplicity might be ~ 30 — 50%, 
based on the positions of brown dwarfs in colour-magnitude 
diagrams (Pinfield et al., 2003), and the statistics of radial 
velocity variations (Maxted and Jeffries, 2005). However, 
Joergens (2006a) has examined the radial velocities of ten 
BDs and very low-mass H-burning stars, and finds no bi- 
naries with separations less than ^0.1 AU, and only two 
objects with variability on timescales greater than 100 days 
that might indicate companions at greater than <~ 0.2 AU. 
Thus, the peak in the separation distribution for BD-BD bi- 
naries is likely to be at ~ 1 to 4 AU. In contrast, the dis- 
tribution of separations in binaries with Sun-like primaries 
peaks at ~ 30 AU {Duquennoy & Mayor, 1991). The mass 
ratio distribution also seems to be dependent on primary 
mass, with BD-BD binaries having a distribution which 
peaks towards equal masses (q PEAK ~ 1), whilst binaries 
with Sun-like primaries have g PEAK ~ 0.3. The implica- 
tion is that, as the primary mass decreases, (i) the multiplic- 
ity decreases, (ii) the distribution of semi-major axes shifts 
to smaller separations and becomes narrower (logarithmi- 
cally), and (iii) the distribution of mass ratios shifts towards 
unity - with these trends all continuing across the divide 
between H-burning stars and brown dwarfs. 

Exotica. Finally, we note that some brown dwarfs are 
components in more complex systems. There are at least 
six systems currently known in which an H-burning star is 
orbited at large distance (> 50 AU) by a BD-BD binary. 
Indeed, preliminary results suggest that wide brown-dwarf 
companions to H-burning stars are 2 to 3 times more likely 
to be in a close BD-BD binary than are field brown dwarfs 
(Burgasser et al, 2005). Cases also exist in which a close 
binary with H-burning components is orbited by a wide 
brown dwarf, while Bouy et al. (2005) have recently re- 
ported the discovery of what is likely to be a triple brown 
dwarf system. However, the statistics of these exotic sys- 
tems are currently too small to interpret quantitatively. 

2.3. Discs, accretion and outflows 

Young brown dwarfs are observed to have infrared ex- 
cesses indicative of circumstellar discs, like young H- 
burning stars (Muench et al, 2001; Natta & Testi, 2001; 
Jayawardhana et al, 2003; Mohanty et al, 2004). Brown- 
dwarf disc lifetimes are estimated to be 3 to 10 Myr, again 
like H-burning stars. From their Ha emission-line pro- 
files, there is evidence for ongoing magnetospheric ac- 
cretion onto brown dwarfs (Scholz and Eisloffel, 2004), 
and the inferred accretion rates form a continuous distribu- 
tion with those for H-burning stars, fitted approximately by 
M ~ 1O- 8 M yr" 1 (M/M Q ) 2 (Muzerolle et al, 2003, 



2005). To sustain these estimated accretion rates, brown 
dwarfs only require rather low-mass discs <~ 1O~ 4 M . 
Finally, the spectra of brown dwarfs also show forbidden 
emission lines suggestive of outflows like those from H- 
burning stars (Fernandez and Comeron, 2001; Natta et al, 
2004), and recently an outflow from a brown dwarf has been 
resolved spatially (Whelan et al., 2005). Thus, in the details 
of their circumstellar discs, accretion rates and outflows, 
young brown dwarfs appear to mimic H-burning stars very 
closely, and to differ significantly only in scale. 

2.4. Rotation and X-rays 

The rotational properties of brown dwarfs also appear to 
connect smoothly with those of very low-mass H-burning 
stars (e.g. Joergens et al, 2003). There is a decrease in the 
amplitude of periodic photometric variations with decreas- 
ing mass, presumably because the decreasing surface tem- 
perature leads to weaker coupling between the gas and the 
magnetic field and hence smaller spots. As with H-burning 
stars, brown dwarfs show evidence for braking by their ac- 
cretion discs; those with strong accretion signatures (broad 
Ha emission) are exclusively slow rotators. 

X-ray emission is detected from M-type brown dwarfs, 
with properties very similar to the X-ray emission from very 
low-mass H-burning stars (Preibisch et al, 2005), but none 
is detected from L-type brown dwarfs, which are too cool 
to support surface magnetic activity. 

2.5. Synopsis 

Given the continuity of statistical properties between 
brown dwarfs and H-burning stars, it is probably unhelp- 
ful to distinguish the formation of brown dwarfs from the 
formation of stars, and in the rest of this review we will 
only use the H-burning limit at ~ 0.075M Q as a reference 
point in the range of stellar masses. The D-burning limit at 
<~ 0.012 M falls in the same category. We will then de- 
fine a star as any object forming on a dynamical timescale, 
by gravitational instability. With this definition, there may 
well be a small overlap between the mass range of stars and 
that of planets. Given that in the immediate future we are 
unlikely to know too much more than the masses and radii 
of the lowest-mass objects, and certainly not their internal 
composition, we will simply have to accept that there is a 
grey area in the range 0.001 to 0.01 M which may harbour 
both stars and planets, and possibly even hybrids. 

It follows that understanding how brown dwarfs form is 
important, not just for its own sake, but because it is a key 
element in understanding why most stars have masses in the 
range 0.01 M to 100 M - and hence why there are lots 
of hospitable stars like the Sun with long-lived habitable 
zones, and enough heavy elements to produce rocky planets 
and life. The high-mass cut-off is probably due to the fact 
that radiation pressure makes it hard to form the highest- 
mass stars. The low-mass cut-off is probably due to the 
opacity limit. By studying brown dwarf formation we seek 
to confirm and quantify the low-mass cut-off. 
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3. STAR FORMATION THERMODYNAMICS 

In this section we review the basic thermodynamics of 
gravitational collapse and fragmentation. The first three 
subsections deal - respectively - with 3-D collapse and hi- 
erarchical fragmentation (§3.1); 2-D one-shot fragmenta- 
tion of a shock compressed layer (§3.2); and fragmentation 
of a disc (§3.3). We describe and contrast these different 
environments, and in each case estimate the minimum mass 
for star formation. In §3.3 we conclude that brown dwarf 
formation by fragmentation of discs around Sun-like stars 
is more likely in the cooler outer parts (R > 100 AU), and 
this may explain why brown dwarf companions to Sun-like 
stars almost always have large separations. In §3.4 we ex- 
plain why impulsive compression does not promote cooling 
of an optically thick fragment. In §3.5 we suggest that close 
BD-BD binaries may be formed by secondary fragmenta- 
tion promoted by the softening of the equation of state when 
H 2 dissociates, and enhanced cooling due to the opacity 
gap. In §3.6 we speculate that close BD-BD binaries may 
form in the outer parts of discs around Sun-like stars. 

We caution that analytic estimates cannot capture all the 
non-linear effects which are likely to occur, and are proba- 
bly important, in a process as chaotic as star formation; they 
are therefore only indicative. A full understanding of any 
mode of star formation requires detailed simulations with 
all potentially influential physical effects included. How- 
ever, as long as converged, robust simulations with all this 
physics properly included remain beyond the compass of 
current supercomputers, analytic estimates provide useful 
insights into the trends to be expected. 

In §3.1 to §3.3, we will be mainly concerned with con- 
temporary star formation in the disc of the Milky Way, and 
therefore with molecular hydrogen at temperatures T < 
100 K where the rotational levels are not strongly excited. 
In this regime the adiabatic exponent is 7 ~ 5/3 and 
the isothermal sound speed is a ~ 0.06 kms -1 (T/K) 1 / 2 . 
With these assumptions, our estimates will not apply to the 
hot gas (T > 10 3 K) where the equation of state is soft- 
ened by effects like H 2 dissociation (e.g. the inner regions 
of discs). We will also assume that the metallicity is ap- 
proximately solar, and that the Rosseland- and Planck-mean 
opacities due to dust are (to order of magnitude) the same, 
R K (T) ~ k p (T) ~ K^T/Kf, with k ± = 10" 3 cm 2 g" 1 
and emissivity index (3 — 2 in the far-infrared and submil- 
limetre. With this assumption, our estimates will again not 
apply to the hot gas (T > 10 3 K) where dust sublimates 
and the opacity falls abruptly with increasing temperature 
(before picking up at even higher temperatures, due to the 
H ion). In §3.4 to §3.6, we relax these assumptions. 

3.1. 3-D collapse and hierarchical fragmentation 

In a uniform 3-D medium, an approximately spherical 
fragment of mass M F3 will only condense out if it is suffi- 
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(Subscript F is for fragment, J is for Jeans, and 3 is for 3-D.) 
The timescale on which the fragment condenses out is 
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The molecular-cloud gas from which stars are forming 
today in the Milky Way is expected to be approximately 
isothermal, with T ~ 10 K, as long as it can radiate effi- 
ciently via molecular lines and dust continuum. Therefore 
it has been argued, following Hoyle (1953), that star forma- 
tion proceeds in molecular clouds by a process of hierarchi- 
cal fragmentation in which an initially massive low-density 
cloud (destined to form a proto-cluster of stars) satisfies 
condition (1) and starts to contract. Once its density has in- 
creased by a factor f 2 , M T3 is reduced by a factor and 
hence parts of the cloud can condense out independently, 
thereby breaking the cloud up into < / sub-clouds. More- 
over, as long as the gas remains approximately isothermal, 
the process can repeat itself recursively, breaking the cloud 
up into ever smaller 'sub-sub.. .sub-clouds'. 

The process ends when the smallest sub-sub. ..sub-clouds 
are so optically thick, and/or collapsing so fast, that the 
PdV work being done on them cannot be radiated away 
fast enough and they heat up. This process is presumed to 
determine the minimum mass for star formation (e.g. Rees, 
1976; Low and Lynden-Bell, 1976), and is usually referred 
to as The Opacity Limit, but see Masunaga and Inutsuka 
(1999) for a more accurate discussion. To estimate M MIN3 
we first formulate the PdV heating rate for a spherical frag- 
ment, neglecting the background radiation field, 
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where in putting dR F3 /dt <~ — (GM F3 /R F3 y/ 2 we are as- 
suming the collapse is dynamical. By comparison, the max- 
imum radiative luminosity of a spherical fragment is 
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where the optical depths are given by r R (T) : 
3M F3Kl (T/K) 2 /4nRl 3 . 

If we follow Rees (1976) and assume that the fragment is 
marginally optically thick, we can put (f R (T) + f p _1 (T)) ~ 
2, and the requirement that C>TL then reduces to 



Pf3 rS Pc3 



.29 



2 2 3 5 5 6 G 3 M 2 3 c 12 h 18 



1/7 



(6) 



4 



Conditions (2) and (6) require p J3 < p C3 and hence 
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Herem PL = (hc/G) 1 ^ 2 = 5.5 x 10~ 5 g is the Planck mass, 
and so M MIN3 is essentially the Chandrasekhar mass times 
a factor (a/c) 1 / 2 ~ 10~ 3 . We also note the relatively weak 
dependence of M MIN3 on T (oc T 1 / 4 ) and the relatively 
strong dependence on m (oc m~ 9 / 4 ). For contemporary 
local star formation we substitute m ~ 4.0 x 10~ 24 g and 
a ~ 1.8 x 10 4 cms" 1 to obtain M MIN3 - 0.004 M . 

In general, the limiting fragment will not necessarily be 
marginally optically thick, but it is trivial to substitute for 
f R (T) and f p (T), and it turns out that, for contemporary 
star formation, the value of M MIN3 is unchanged. This is 
because - coincidentally - the limiting fragment in contem- 
porary star formation is marginally optically thick. 

There are, however, some serious problems with 3-D hi- 
erarchical fragmentation. There is no conclusive evidence 
that it operates in nature, and nor does it seem to occur in 
numerical simulations of star formation. This is because a 
proto-fragment inevitably condenses out more slowly than 
the larger structure of which it is part, by virtue of the fact 
that it is, at every stage, less Jeans unstable than that larger 
structure (see Eqn. 3). Therefore the proto-fragment is 
very likely to be merged with other nearby fragments be- 
fore its condensation becomes nonlinear. In addition, even 
if a proto-fragment starts off with mass ~ M J3 , it will sub- 
sequently increase its mass by a large factor through accre- 
tion, before its condensation becomes nonlinear. Finally, 
individual fragments will be back-warmed by the ambient 
radiation field from other cooling fragments, which in prin- 
ciple fill a significant fraction of the celestial sphere, and 
again this will tend to increase M MIN3 . 

3.2. 2-D one-shot fragmentation of a shocked layer 

In fact, 3-D hierarchical fragmentation may be an inap- 
propriate paradigm for star formation in molecular clouds. 
There is growing evidence that, once a molecular cloud is 
assembled, star formation proceeds very rapidly, essentially 
'in a crossing time' (Elmegreen, 2000). In this scenario, star 
formation occurs in molecular clouds only where two or 
more turbulent flows of sufficient density collide with suffi- 
cient ram pressure to produce a shock-compressed layer or 
filament which then fragments to produce prestellar cores; 
in cases where flows converge simultaneously from several 
directions, isolated cores may even form (see Padoan and 
Nordlund, 2002, 2004; and §4). 

A basic model of this scenario can be constructed with 
relatively few parameters by considering two flows having 
pre-shock density p, colliding at relative speed v, to produce 
a shock-compressed layer. If the effective isothermal sound 
speed in the resulting layer is a, the density is ~ p(v/a) 2 . 
The layer is initially contained by the ram pressure of the in- 
flowing gas, and until it fragments it has a rather flat density 
profile. It fragments at time t F2 , whilst it is still accumulat- 



ing, and the fastest growing fragments have mass M F2 , ra- 
dius R F2 (in the plane of the layer) and half-thickness Z F2 
(perpendicular to the plane of the layer) given by 
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(see Whitworth et al, 1994a, 1994b). We note (a) that 
this mode of fragmentation is '2-D' because the motions 
which assemble a fragment out of the shock-compressed 
layer are initially largely in the plane of the layer; (b) 
that it is 'one-shot' in the sense of not being hierarchi- 
cal; (c) that the fragments are initially flattened objects 
(R F2 /Z F2 ~ v/2a > 1); (d) that M F2 is not simply the 
standard 3-D Jeans mass (M J3 , Eqn. 1) evaluated at the 
post-shock density and velocity dispersion - it is larger by 
a factor ~ (v/a) 1 / 2 ; and (e) that in reality the colliding 
flows will contain density substructure which acts as seeds 
for fragmentation and gives rise to a range of M F2 values. 

From Eqn. (9), we see that the fragment mass decreases 
monotonically with the mass-flux in the colliding flows, pv. 
As in hierarchical fragmentation, a fragment in a shock- 
compressed layer will only condense out if it is able to re- 
main approximately isothermal by radiating efficiently. The 
PdV heating rate for a flattened fragment in a layer is 
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The radiative cooling rate of the fragment is 

8ir 5 fh 4 a n /15c 2 h 3 Gpv 



2 7 ri? 2 2 a SB T 4 



(f R (T)+f p -i(T)) (r R (T)+f p -i(T)) 



(13) 

where the optical depths are now given by r R (T) = 
f p (T) = (2apv /irG) 1 ^ 2 k 1 (T/K) 2 . The requirement that 
C > TL then reduces to a limit on the mass flux in the 
colliding flows, 
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If we assume the fragment is marginally optically thick, and 
set (f R (T) + V 1 ^)) - % we obtain 
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in: 
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and for contemporary local star formation, this gives 
M MIN2 ~ 0.001 M Q . Once again, if we treat the com- 
pletely general case by including the optical-depth terms, 
we obtain essentially the same value for M MIN2 , because - 
purely coincidentally - the limiting mass for contemporary 
local star formation is again marginally optically thick. 

Although one-shot 2-D fragmentation of a shock- 
compressed layer and 3-D hierarchical fragmentation give 
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essentially the same expression for the minimum mass 



3 ), they give very different expressions for 



the Jeans mass (M J2 ^ Af, 3 ), and there are other impor- 
tant differences. In particular, 2-D one-shot fragmentation 
bypasses all the problems associated with 3-D hierarchi- 
cal fragmentation. There is no backwarming because there 
are no other local fragments filling the part of the celestial 
sphere into which a fragment radiates (i.e. perpendicular 
to the layer). More importantly, there is little likelihood of 
fragments merging, since fragments with M F2 ~ M J2 con- 
dense out faster than any other structures in a layer (whereas 
fragments in a 3-D medium with M F3 ~ M 13 condense out 
slower than all larger structures in that medium; e.g. Lar- 
son, 1985). Finally, because condensation in a layer is so 
fast, growth by accretion is limited. 

Boyd and Whitworth (2004) have analysed the radiative 
cooling of a fragmenting layer, taking into account not only 
the PdV heating of a condensing fragment, but also the on- 
going accretion (as matter continues to flow into the layer) 
and the energy dissipated in the accretion shock. They find 
that for T ~ 10 K, the minimum mass is 0.0027 M . This 
fragment starts with mass ~ 0.0011 M , but continues to 
grow by accretion as it condenses out. 

3.3. Fragmentation of a circumstellar disc 

Another scenario which may be more relevant to con- 
temporary star formation than 3-D hierarchical fragmenta- 
tion is fragmentation of a circumstellar disc. There are three 
critical issues here, (i) Under what circumstances does a 
disc fragment gravitationally? (ii) Can the resulting frag- 
ments cool fast enough to condense out? And (iii) can the 
resulting fragments lose angular momentum fast enough to 
condense out? The last two issues are critical because, if 
a fragment cannot cool and lose angular momentum fast 
enough, it is likely to bounce and be sheared apart. For 
simplicity we consider an equilibrium disc. 

(i) The condition for an isolated disc to fragment gravi- 
tationally is that the surface density, E, be sufficiently large, 



E > XL ~ 



a e 
7rG' 



(16) 



where a is the local sound speed and e is the local epicyclic 
frequency (Toomre, 1964). 

Condition (16) is also the condition for spiral modes to 
develop in the disc, and these will have the effect of redis- 
tributing angular momentum. As a result, the inner parts 
of the disc may simply accrete onto the central primary 
star and the outer parts may disperse without fragmenting 
(Laughlin and Bodenheimer, 1994; Nelson et al, 1998). 
Thus, if fragmentation is to occur, it must occur on a dy- 
namical timescale. 

(ii) The condition for a fragment to cool fast enough 
to condense out is therefore that the fragment can radiate 
away, on a dynamical timescale, the thermal energy being 
delivered by compression. Gammie (2001) has shown that 
for a Keplerian disc, this condition can be written as a con- 



straint on the cooling time, 

t 



< 1 

COOL r^j Q 1 



(17) 



where il is the local orbital angular speed. 

We assume that, as a disc forms, E increases sufficiently 
slowly that it does not greatly exceed E T when the disc be- 
comes unstable. It then follows that the radius, growth time 
and mass of the fastest growing fragment are 



The compressional heating rate for a fragment is thus 

H-P dV _ 3MfdQ2 _ 3 " 5 - 
dt 2t™ 2G ' 



(18) 
(19) 
(20) 

(21) 



and the radiative cooling rate is 



2nR 2 a SB T 4 

r ^ fd ^ B 



47r 6 m 4 a 10 /15c 2 fe 3 e 2 
(r R (T)+f p -i(T)) " (f R (T) + f p -i(T)) 



, (22) 



where now r R (T) ~ r p (T) ~ a€K(T)/irG. Consequently 
the requirement that C>H reduces to 



e 2 47r 6 Gm 4 /45c 2 /i 3 
^~ (f R (T)+f p -i(T)) 



(23) 



To illustrate the discussion we consider the specific case 
of a Keplerian disc around a Sun-like star, with 



M, 



AU 



T(D) * 300K(^ 



(24) 



(25) 



and hence a(D) ~ lkms" 1 (L±/L Q ) 1/8 (D/AU) _1/4 , 
where D is distance from the Sun-like star. The fastest 
growing fragment then has mass 



M FD *3xlO-«M 0VM ^ u ^ 



M*V V2 {LA' V ° f_D_\ 



3/8 



3/4 



and Condition (23) is only satisfied for 



D > 90 AU 



M FD > 0.003 M 





( — 






/AM 


-5/28 







-1/7 



L 



15/56 



(26) 



(27) 



(28) 



(iii) Angular momentum is removed from a condens- 
ing fragment by gravitational torques, in the same way that 
angular momentum is redistributed in the disc as a whole 
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when E > S T . This is a non-linear and stochastic process, 
and there is no analytic estimate of the rate at which it oc- 
curs. Therefore the condition for a fragment to lose angular 
momentum fast enough to condense out has to be evalu- 
ated by numerical simulation. The chapter by Durisen et al. 
deals with this problem. 

3.4. Non-linear thermodynamics 

Impulsive triggers which produce rapid compression 
will always help to amplify self-gravity, because the freefall 
time varies as t FF oc p~ x / 2 . However, there is only a very 
restricted temperature range within which rapid compres- 
sion will help an optically thick proto-fragment to cool 
more rapidly (and thereby avoid the likelihood of its bounc- 
ing and being sheared apart). Suppose that the Rosseland- 
mean opacity is given by k b oc p a T 13 , and that the gas has 
a ratio of specific heats 7. 

Then if a fragment condenses out quasistatically, its 
cooling time varies as i COOL oc p a +( 1 +' 3 )/ 3 , so i COOL 
only decreases as fast as t FF , with increasing density, if 
(i < — (6a + 5)/2. This condition is only likely to be sat- 
isfied in the temperature range where refractory dust subli- 
mates (1500 to 3000 K). 

If instead the fragment is compressed impulsively - and 
therefore adiabatically - its cooling time varies as i COOL oc 
pa+4/3+(/3-3)(7-i) > so now t COOL only decreases as fast as 
t FF , with increasing density, if/3 < 3 — (6a + 1 1) /6 (7 — 1 ) . 
For T < 100 K, we have 7 ~ 5/3 and so enhanced cool- 
ing requires (3 < —1/4, which is unlikely. For 100 K < 
T < 1000 K, 7 ~ 7/5 and enhanced cooling requires 
P % -19/12, which is even less likely. For 1,000K < 
T < 3, 000 K, H 2 dissociation gives 7 ~ 1.1 and enhanced 
cooling requires (3 < —15, which may occur during sub- 
limation, but then necessarily only over a small tempera- 
ture range. For 3,000K < T < 10, 000 K, 7 - 5/3 and 
H opacity gives a <~ 1/2 and [3 ~ 4, so t COOL oc p 1 ! 2 . 
At even higher temperatures where Kramers opacity dom- 
inates t COOL oc p~ 2 , and where electron scattering opacity 
dominates t COOL oc p~ 2 ^ 3 , but by this stage a fragment is 
very opaque and very strongly bound. 

3.5. Forming close BD-BD binaries 

It is possible that there is a secondary fragmentation 
regime at T ~ 2000 K, due to the softening of the equa- 
tion of state caused by H 2 dissociation, and the enhanced 
cooling which occurs in the opacity gap between dust sub- 
limation and H opacity. A spherical cloud of mass M in 
hydrostatic equilibrium at this temperature has radius 

R ~ 4AU (oft) ' (29) 

If the cloud fragments into a binary and virialises, the binary 
should have a separation of the same order. Eqn. (29) is 
actually a good mean fit to the values plotted on Fig. 6 of the 
chapter by Burgasser et al, suggesting that BD-BD binaries 



may be produced by secondary fragmentation facilitated by 
H 2 dissociation and/or the opacity gap. 

3.6. Forming close BD-BD binaries in discs 

Burgasser et al. (2005) have noted that - modulo the 
small-number statistics involved - a brown dwarf in a wide 
orbit (> 200 AU) about a Sun-like star is apparently more 
likely to be in a close BD-BD binary system (< 20 AU) 
than a brown dwarf in the field. If this trend is confirmed, it 
suggests that BD-BD binaries are formed in discs, and then 
may be ejected. For a close BD-BD binary in a wide orbit 
around a Sun-like star, the internal binding energy of the 
BD-BD binary is typically comparable with or larger than 
the binding energy of the BD-BD binary to the Sun-like 
star, and therefore some BD-BD binaries should be able to 
survive ejection. 

Moreover, we know, from the Toomre criterion, that for 
fragments condensing out of discs, losing angular momen- 
tum is a critical issue. The smaller a condensation becomes, 
the slower the rate at which angular momentum can be lost 
by gravitational torques. Therefore such a condensation 
may be strongly disposed to binary fragmentation when 
its thermal support is weakened at T ~ 2000 K. Since 
only a subset of brown dwarfs is in close BD-BD binaries, 
this suggestion does not require that all brown dwarfs are 
formed in discs. 

4. FORMING VERY LOW-MASS PRESTELLAR 
CORES IN TURBULENT CLOUDS 

A number of surveys (Testi and Sargent, 1998; Peng et 
al, 1998; Motte et al, 1998; Motte et al, 2001; Motte and 
Andre, 2001; Johnstone et al, 2000; Sandell and Knee, 
2001; Onishi et al, 2002; Reid and Wilson, 2005) have 
noted the close similarity between the core mass function 
(CMF) and the stellar initial mass function (IMF). Cores 
with brown dwarf masses have been observed in several 
of these surveys, but since these low-mass cores are usu- 
ally below the completeness limit, their statistics are un- 
reliable. However, the similarity between the CMF and the 
IMF suggests that each core gives rise to a similar low num- 
ber of stars, whose final masses are heavily influenced by 
the mass of the core. Indeed, protostars are generally ob- 
served inside relatively well defined envelopes, as implied 
by the definition of Class and Class I objects (Andre et 
al, 1993, 2000; Tachihara et al, 2002). This is not to say 
that the final masses of stars are the same as the masses of 
the prestellar cores in which they form. Protostars are asso- 
ciated with outflows, which can remove a significant frac- 
tion of the mass of the protostellar envelope, and cores may 
split into more than one star (see also Hosking and Whit- 
worth, 2004; Stamatellos et al, 2005). However, low mass 
cores appear typically to harbour only one or two protostars 
(e.g. Tachihara et al, 2002), and the overall binary statistics 
are inconsistent with a larger number of stars being formed 
from each core (Goodwin and Kroupa, 2005). 
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4.1. Core structure and low mass core formation 

Observed prestellar cores generally have Bonnor-Ebert 
(BE) like density profiles, with relatively sharp outer edges 
(Bacmann et al., 2000; Motte and Andre, 2001; Kirk et al., 
2005). Starless cores have flat BE-like profiles, while cores 
with detected Class or Class I protostars have more cen- 
trally peaked density profiles, presumably as a result of their 
deeper potential wells. Since well defined cores with sharp 
edges define finite mass reservoirs for the stars that form 
inside them, the core observations reinforce the notion that 
the masses of stars are, at least statistically, strongly influ- 
enced by the masses of the cores within which they form 
(Padoan and Nordlund, 2002, 2004). 

It is therefore important to consider whether it is possible 
to form brown dwarfs directly, from correspondingly low- 
mass cores. The standard argument against brown dwarfs 
being formed directly is that the density needs to be very 
high for a fragment with brown dwarf mass to collapse 
gravitationally. If one considers low-amplitude perturba- 
tions on top of a typical mean density of 10 4 to 10 5 cm~ 3 , 
he will conclude that only fluctuations of order the Jeans 
mass (1 to 3M ; Eqn. 1) are able to collapse in a typical 
molecular cloud. However, many cores actually have very 
high density contrast relative to their surroundings, and so 
it is inappropriate to use the Jeans mass at the mean density 
as an estimate of the resulting protostellar mass. 

It is also important to realise that cores having suffi- 
ciently high density to form brown dwarfs directly need not 
- indeed must not - be very common, and therefore we can 
appeal to exceptional circumstances to generate them. As 
a measure of how exceptional the circumstances must be, 
consider the IMF (in the form dN / d In AO- At high masses, 
the IMF falls off with the Salpeter (1955) exponent -1.35, 
and if this power-law IMF continued unbroken down to the 
brown-dwarf regime, I0 4x1 - 35 ~ 250, 000 brown dwarfs 
with M ~ 0.01 M would be formed for every massive star 
with M ~ 100 M . Instead, the IMF peaks near ~ 0.3 M 
and falls off at low masses; adopting the Chabrier (2003) 
IMF for the Milky Way, we find that the number of brown 
dwarfs formed with M ~ 0.01 M is actually about equal 
to the number of massive stars formed with M ~ 100 M , 
as predicted theoretically by Zinnecker (1984). 

In order for a high-density core to form, its mass must be 
concentrated into a small volume. If u is the fluid velocity, 
accumulation of mass is measured by the quantity — V • u = 
D\np/Dt (the comoving time derivative of log density), 
and so individual cores form at local maxima of — V • u. 

The flow towards a convergence point is generally super- 
sonic, and so stand-off shocks develop, separating the un- 
shocked upstream gas from the shocked and nearly stagnant 
downstream gas. The stand-off shocks correspond to den- 
sity jumps ~ M 2 if the shocked gas is dominated by gas 
pressure; and ~ A4 A (the Alfvenic Mach number) if the 
shocked gas is dominated by magnetic pressure. As mat- 
ter accumulates around the convergence point, a deepening 
gravitational potential well develops, and the accumulated 



and stagnant gas forms a growing Bonnor-Ebert like core, 
stabilised by the external ram pressure, /o|u| 2 . 

If the density at the centre of the core increases to more 
than about 14 times the boundary density (just inside the 
shock surface), the core becomes gravitationally unstable 
and collapses. However, even before it becomes unstable, 
such a core has - by virtue of the shock jump at its bound- 
ary - a total density contrast (between its centre and the 
surrounding inflowing unshocked gas) that can greatly ex- 
ceed a factor of 14. Moreover, if the inflowing gas runs out 
before the core has become sufficiently massive and dense 
to collapse, then the core will expand and disperse due to 
the decrease in ram pressure at its boundary. 

In the last column of Table 3 in Bacmann et al. (2000), at 
least 8 of the 9 cores listed have density contrasts exceed- 
ing 14, but no strong indication of collapse. Bacmann et 
al. conclude that the structures are best fitted with Bonnor- 
Ebert spheres, but since they cannot explain how the cores 
can be stable at these high density contrasts, they go on to 
discuss other density profiles and models, which do not fit 
the data as well as BE-spheres. Similar core profiles have 
been observed in many other surveys. There is thus direct 
observational evidence for the creation of prestellar cores 
with large density contrasts between their centres and the 
surrounding medium, as expected when cores are produced 
as stagnant structures by supersonic flows. Indeed, cores 
formed in numerical simulations of molecular cloud turbu- 
lence have internal velocity dispersions and rotation veloci- 
ties that are entirely consistent with those of observed cores 
(cf. Figs. 8 and 9 in Nordlund and Padoan, 2003). 

With this scenario in mind, we now investigate semi- 
quantitatively the range of densities that is possible, and the 
circumstances under which the densities attained would be 
high enough to lead to the collapse of brown dwarf mass 
cores. Any particular converging flow may be broadly char- 
acterised by three parameters: the up-stream density p , the 
Mach number M of the up-stream flow, and the degree of 
focussing of the up-stream flow towards a 3-D convergence 
point. As a measure of the latter we take the ratio / of the 
surface area that encloses the up-stream flow to the surface 
area of the stand-off shock that surrounds the central BE- 
like structure. Since the up-stream flow is supersonic, its 
motion is essentially inertial (i.e. constant velocity), until 
it encounters the stand-off shock. Thus, by mass conser- 
vation, its density increases from p to fp . At the stand- 
off shock the density increases by an additional factor M 2 
for hydro-shocks (or Ai A for MHD-shocks). Finally, from 
just inside the shock to the centre of the BE-like core the 
density increases by a further factor < 14 (~ 14, in the 
marginally stable case). The total density increase from the 
up-stream source to the core centre is ~ 14/.M 2 . This fac- 
tor has the right dependence to account for the formation of 
brown dwarf mass cores, in that it can be arbitrarily large, 
albeit with decreasing probability. 

As an example, consider the formation of a core with 
temperature ~ 10 K and central density ~ 10 7 cm~ 3 , suffi- 
ciently high to form a star with M ~ 0.1 M . If the average 
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Fig. 1. — Density- and velocity-profiles, through a core forming 
from a spherically symmetric convergent flow. The thin line rep- 
resents an early stage when the ~ 0.015 M core approximates to 
a stable and stagnant BE-like configuration; the boundary shock 
is at ~ 300 AU. The thick line represents a later stage after the 
core has become unstable and started to collapse; the boundary 
shock is now at ~ 350 AU, because the core is more massive 
(~ 0.025 M ). Outside the shock is the convergent inertial flow. 
Note the very large density contrast, even between the gas at the 
edge of the core and the ambient medium at ~ 1500 AU. 

cloud density is <~ 10 4 cm~ 3 and the Mach number is 5, / 
only needs to be ~ 3; i.e., the up-stream flow needs to be 
focussed onto a three times smaller stand-off shock surface. 
This must be rather common, and indeed at M ~ 0.1 M 
the IMF has hardly dropped below its maximum. 

Repeating this estimate for a core with M ~ 0.01 M 
one finds that the flow needs to be much better focussed. 
The stand-off shock surface now needs to be ~ 300 times 
smaller than the up-stream source area, and so the linear 
size of the core needs to be ~ 17 times smaller than the 
scale of the up-stream flow. This does require a rather ex- 
ceptional focussing of the up-stream flow, but then brown 
dwarfs with M ~ 0.01 M actually are very rare. The frac- 
tion of brown dwarfs that form in this manner can only be 
ascertained quantitatively by performing very high resolu- 
tion numerical simulations. 

Fig. 1 illustrates an example intermediate between these 
two extremes, i.e. a core with final mass ~ 0.03 M . 

4.2. Envelope break-out and continued accretion 

The sharp density transitions surrounding accreting 
cores produce more well defined cores than would exist 
in a subsonic medium. Nevertheless there will be some 
continued accretion after a collapsed object is formed, until 
the upstream supply of gas is exhausted. Except towards 
the end of the process, such cores correspond to Class 
objects, which by definition have not yet accreted half their 
final mass (Andre et al, 2000). 

Class I and II objects are observed to accrete from low 
mass envelopes (Motte and Andre, 2001), and have much 
reduced accretion rates, < 10~ 6 M yr _1 . Once the natal 



envelope is consumed there is no longer a strong coupling 
between the collapsed object (moving under the N-body in- 
fluence of neighbouring objects and gas) and the surround- 
ing medium, and one expects the collapsed object to pick 
up speed relative to the surrounding gas, and to then accrete 
in a manner similar to Bondi-Hoyle accretion (Padoan et 
al., 2005; Krumholz et al, 2005b). In a turbulent medium 
the problem is more complicated {Krumholz et al, 2005a) , 
but the scaling is more or less as for Bondi-Hoyle accretion, 
with a prefactor $ ~ 1 to 5 that accounts for the compli- 
cations. Quantitative estimates indicate that the accretion 
after break-out from the natal core is insignificant. 

4.3. The viability of turbulent fragmentation 

Although turbulent fragmentation generates a mass func- 
tion for prestellar cores, matching broadly the observed stel- 
lar IMF, there are two caveats which should be born in mind 
when considering the formation of brown dwarfs. 

First, in turbulent fragmentation (Padoan and Nordlund, 
2002), brown dwarfs are produced by a small subset of low- 
mass cores, i.e. those which are sufficiently dense to col- 
lapse gravitationally. At its low-mass end, the overall (i.e. 
fully sampled) CMF should be dominated by a much larger 
number of transient (i.e. non-prestellar) cores. Moreover, 
many of these transient cores will be only slightly less dense 
than the ones that spawn brown dwarfs. They will therefore 
be detectable in submillimetre continuum observations, and 
they should also be somewhat longer lived. Thus, even al- 
lowing for selection effects, the observed CMF should fall 
much less steeply with decreasing mass than the stellar IMF. 
Recent estimates of the CMF (Nutter & Ward-Thompson, in 
prep.), which probe to lower masses by using longer expo- 
sures, actually suggest the opposite, but completeness re- 
mains a concern (e.g. Kirk et al., 2006). 

Second, turbulent fragmentation predicts a ratio of 
brown dwarfs to H-burning stars, 1Z, which is exponen- 
tially sensitive to the Alfvenic Mach number on the largest 
scales (Ai A ) and to the mean cloud density (n). Regions 
with smaller values of M A and/or n should generate stel- 
lar populations with significantly fewer brown dwarfs, and 
regions with larger values should generate stellar popula- 
tions with significantly more brown dwarfs. However, in 
nature 1Z appears to vary little over a wide range of local 
star-forming environments (see chapter by Luhman et al). 

This problem can be overcome if nature selects a narrow 
range of M A and n. For example, Whitworth (2005) has 
noted that contemporary star formation may only proceed 
rapidly if the gas couples thermally to the dust (so that it can 
avail itself of broadband - as distinct from molecular-line - 
cooling). This requires that the ram pressure in shocks pro- 
ducing prestellar cores exceeds a critical value, which con- 
verts into the constraint nTM 2 A > P CKIT ~ 10 5 cm~ 3 K. 
This constraint may help to select the combinations of A4 A 
and n which reproduce the observed ratio of brown dwarfs 
to H-burning stars. 
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5. COLLAPSE AND FRAGMENTATION OF LARGE 
PRESTELLAR CORES 

Whilst a very low-mass prestellar core (< 0.1 M ) must 
collapse to form either a single brown dwarf or a multiple 
brown dwarf system, larger prestellar cores (> 1 M ) are 
expected to form clusters of stars having a range of masses. 
We can identify five mechanisms which may play a role 
in determining the final stellar masses, (i) During the ap- 
proximately isothermal initial collapse phase, as the pres- 
sure becomes increasingly unimportant and the collapse 
approaches freefall, self-gravity will amplify any existing 
density substructure, (ii) Then, when the density reaches 
n H2 ~ 10 11 cm~ 3 , the gas becomes optically thick and 
switches rather suddenly from approximate isothermality 
to approximate adiabaticity. At this juncture, a network of 
shock waves develops to slow the collapse down, and non- 
linear interactions between these shock waves produce and 
amplify further substructure (sheets, filaments and isolated 
prestellar cores), (iii) Some of these structures will have 
sufficient angular momentum to form discs, and these may 
then fragment due to rotational instability (see §6). Finally, 
once a protostellar embryo (i.e. an object which is suffi- 
ciently well bound to be treated dynamically as a single en- 
tity) has condensed out of a fragment, its subsequent evolu- 
tion and final mass will be determined by (iv) competitive 
accretion and (v) dynamical interaction (see §7). Here we 
concentrate on mechanisms (i) and (ii), since these are the 
ones which distinguish the evolution of a high-mass core 
from the low-mass cores considered in §4. 

High-mass prestellar cores (> 10 M ) invariably display 
non-linear internal density structure, and at the typical den- 
sity (nn 2 ~ 10 5 cm~ 3 ) and temperature (T ~ 10 K) in 
a prestellar core, the Jeans mass is M J3 ~ 0.8 M (Eqn. 
1). Therefore, in the absence of a significant magnetic field, 
they are very likely to fragment during collapse. 

Even the smallest cores are usually far from spherical, 
and have been modelled as being either prolate or oblate 
(Myers et al, 1991; Ryden, 1996; Jones et al, 2001; Good- 
win et al, 2002; Curry, 2002; Myers, 2005) with prolate 
models being favoured statistically, although in reality cores 
are probably triaxial - or even more complicated - in their 
full three dimensional structure (Jones et al, 2001; Good- 
win et al, 2002). Gravity works to enhance anisotropies in 
collapsing objects (Lin et al, 1965) with collapse occurring 
fastest along the shortest axis to form sheets and filaments, 
which then subsequently fragment (Bastien, 1983; Inutsuka 
and Miyama, 1992). Thus, it is unsurprising that hydro- 
dynamical simulations of both oblate and prolate cores are 
prone to fragmentation (e.g. Bastien et al, 1991; Boss, 
1996). 

Prestellar cores also tend to have complex internal veloc- 
ity fields (Larson, 1981; Myers and Benson, 1983; Arquilla 
and Goldsmith, 1985), but since prestellar cores can only be 
observed from a single direction, the interpretation of these 
velocities is difficult. If cores are assumed to be in solid- 
body rotation, the ratio of rotational to gravitational energy 



is typically /3 = TZ/\Sl\ ~ 0.03, with some cases as high as 
P ~ 0.1 (Goodman et al, 1993; Barranco and Goodman, 
1998). However, the observed velocities are more likely to 
be turbulent in nature, i.e. less well ordered than solid-body 
rotation (Myers and Gammie, 1999; Burkert and Boden- 
heimer, 2000). Quite low levels of turbulence (e.g. Good- 
win et al. 2004a), and/or global rotation (e.g. Boss, 1986; 
Bonnell and Bate, 1994a; Hennebelle et al, 2004; Cha and 
Whitworth, 2004) are sufficient to make a collapsing core 
fragment into a small ensemble of protostellar embryos. 

Unfortunately, collapse and fragmentation can only be 
explored by means of numerical simulations, and there is 
a huge and poorly constrained range of admissible initial 
conditions, which makes the extraction of robust theorems 
very hard (see e.g. Hennebelle et al, 2004). Moreover 
almost all simulations to date use barotropic equations of 
state (i.e. P = P(p)). These barotropic equations of state 
are designed to mimic the expected thermal behaviour of 
protostellar gas, but they do not capture the thermal inertia 
effects which become important at the juncture when the 
gas starts to heat up due to adiabatic compression, and it 
appears (Boss et al, 2000) that these thermal inertia effects 
play a critical, deterministic role in gravitational fragmenta- 
tion. Proper treatment of the energy equation and the asso- 
ciated radiative transport (e.g. Whitehouse and Bate, 2006) 
is needed to make these simulations more realistic. 

6. DISC FRAGMENTATION 

We organise our discussion of disc fragmentation under 
the headings (i) isolated, relaxed discs, (ii) unrelaxed discs, 
and (iii) interacting discs. 

6.1. Isolated relaxed discs 

The dynamical fragmentation of isolated relaxed discs is 
discussed in detail in the chapter by Durisen et al. Although 
the emphasis there is on the genesis of planets, the same is- 
sues pertain to the formation of brown dwarfs, viz. Under 
what circumstances do discs become unstable against frag- 
mentation? Is Eqn. (16) a sufficient condition for gravita- 
tional instability (in which case, what is the precise value of 
£ T ) or is it also necessary for E to increase rapidly (in or- 
der to avoid the disc simply being accreted and dispersed by 
torques due to spiral density waves)? Does Eqn. (17) deter- 
mine whether fragments can cool fast enough to condense 
out? What role is played by thermodynamic effects like H 2 
dissociation, dust sublimation, and convection? What are 
the properties of dust in discs, and how well mixed are the 
gas and dust? Does the survival of a fragment depend on its 
ability to lose angular momentum rapidly? Some of these 
issues are only beginning to be investigated. 

Rice et al (2003) have corroborated Eqn. (17) numer- 
ically by performing SPH simulations of disc fragmenta- 
tion with a parameterised cooling law of the form du/dt = 
—uVi/ (3 (where u is the specific internal energy and 51 the 
orbital angular speed). Endemic fragmentation occurs with 
/3 = 3 but not with /3 = 5. However, because this cooling 
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law results in indefinite cooling, whereas there is a limit to 
the supply of rotational energy which can be tapped through 
shock heating, by the time fragmentation occurs the tem- 
peratures have dropped to rather low values, which may be 
hard to realise in nature. 

Rafikov (2005) argues that the opacity of gravitationally 
unstable discs is so high, and the cooling times are therefore 
so long, that fragments can only condense out on a dynam- 
ical timescale in the outer, cooler regions of massive discs. 
(Although the treatment is somewhat different, his conclu- 
sion resonates with the simple analysis we have presented 
in §3.3 and carries the same caveats). 

Johnson and Gammie (2003) point out that the effects 
of opacity may be mitigated in temperature regimes where 
dust sublimates and hence the opacity decreases abruptly 
with increasing temperature. However, dust sublimation ef- 
fects are confined to the temperature range T > 200 K, and 
the most critical effects occur at T > 2000 K, so they are 
probably only relevant to fragmentation close to the central 
star (D < 10 AU). 

Cai et al. (2005) report simulations of disc evolution 
taking radiation transport into account, and conclude that - 
even at low metallicities - the opacity in the hot inner parts 
of a disc is too high to allow fragmentation on a dynamical 
timescale, in agreement with Rafikov (2005). 

In contrast, Boss (2004) presents simulations of disc 
evolution taking radiation transport into account and con- 
cludes that fragments of planetary mass do condense out - 
or at least that gravitationally bound fragments form which 
would subsequently condense out. He argues that the proto- 
fragments in his simulations bypass the effects of high 
opacity by transporting energy convectively. However, it is 
not clear how convection can cool a fragment which is con- 
densing out on a dynamical timescale. The coherent small- 
scale motions he attributes to convection may actually be 
manifestations of a fragment which is unable to cool and 
is bouncing - or will soon bounce - prior to being sheared 
apart. This needs to be investigated further, but the numeri- 
cal complexities are considerable. 

6.2. Unrelaxed discs 

In numerical simulations of the collapse and fragmen- 
tation of intermediate- and high-mass prestellar cores, the 
first single protostars to form usually quickly acquire mas- 
sive circumstellar discs, and secondary protostars then con- 
dense out of these discs. This pattern is common for ro- 
tating cores (e.g. Bonnell, 1994; Bonnell and Bate, 1994a, 
1994b; Turner et al, 1995; Whitworth et al, 1995; Burkert 
and Bodenheimer, 1996; Burkert et al, 1997), for turbulent 
cores (e.g. Bate et al, 2002b, 2003; Goodwin et al, 2004b) 
and for cores which are subjected to a sudden increase in 
external pressure (e.g. Hennebelle et al, 2004). The discs 
thus formed fragment before they have time to relax to an 
equilibrium state. Indeed, the material accreting onto the 
disc is often quite lumpy, and this helps to seed fragmenta- 
tion. Under this circumstance, gravitational fragmentation 



is more likely simply because proto-fragments are launched 
directly into the non-linear regime of gravitational instabil- 
ity, rather than having first to grow through the linear phase. 

Even so, a proto-fragment still has to be able to cool 
and lose angular momentum on a dynamical timescale, if 
it is to condense out, rather than bouncing and then being 
sheared apart. Fragmentation of unrelaxed discs can only be 
studied by means of numerical simulations, and to date no 
simulations of the process have been performed which treat 
properly the energy equation and the associated transport 
of radiation. Since young protostars have relatively high lu- 
minosities, an important consideration will be irradiation of 
the disc by the primary protostar at its centre and any other 
nearby protostars. 

6.3. Interacting discs 

Another way in which disc fragmentation can be trig- 
gered is by an impulsive interaction with another disc, or 
with a naked star. In the dense protocluster environment 
where most stars are presumed to form, such interactions 
must be quite frequent, since many very young protostel- 
lar discs have diameters > 300 AU and the mean separa- 
tion between neighbouring protostars in a typical cluster is 
< 3000 AU. Indeed, ~ 50% of solar-type stars end up in 
binary systems with semi-major axes a < 1000 AU, so the 
notion of a disc evolving in the gravitational field of a sin- 
gle, isolated protostar is probably rather artificial. 

Boffin et al (1998) and Watkins et al (1998a, 1998b) 
have simulated parabolic interactions between two proto- 
stellar discs, and between a single protostellar disc and a 
naked protostar. All possible mutual orientations of spin 
and orbit are sampled, and the gas is assumed to behave 
isothermally, which is probably a reasonable assumption, 
since the discs are large (initial radius 1000 AU) and most 
of the secondary protostars form at large distances (peri- 
astra > 100 AU). The critical parameter turns out to be the 
effective shear viscosity in the disc. If the Shakura-Sunyaev 
parameter is low, a ss ~ 10~ 3 , most of the secondary proto- 
stars have masses in the range 0.001 M to 0.01 M . Con- 
versely, if a ss is larger, a ss ~ 10~ 2 , most of the secondary 
protostars have masses in the range 0.01 M to 0.1 M . 
The formation of low-mass companions is most efficient 
for interactions in which the orbital and spin angular mo- 
menta are aligned; on average 2.4 low-mass companions are 
formed per interaction in this case. If the orbital and spin 
angular momenta are randomly orientated, then on average 
1 .2 companions are formed per interaction. It is important 
that such simulations be repeated, with a proper treatment 
of the energy equation and the associated energy transport, 
to check their fidelity, and to establish whether low-mass 
companions can form at closer periastra. 

7. PREMATURE EJECTION OF PROTOSTELLAR 
EMBRYOS 

The scenario where brown dwarfs form by premature 
ejection is closely linked to - but ultimately independent 
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of - the notion of competitive accretion. In the ejection 
hypothesis, brown dwarfs are simply protostellar embryos 
which get separated from their reservoir of accretable ma- 
terial at an early stage, and so in the context of brown dwarf 
formation it is irrelevant whether there are other stars com- 
peting for this same material. The importance of compet- 
itive accretion for intermediate and high-mass stars is ar- 
gued in the chapter by Bonnell et al. The contentious issue 
is whether protostellar embryos - i.e. the first, very low- 
mass (~ 0.003 M ) high-density star-like objects - exist 
for long enough to do much competing. At one extreme, 
it is argued that a protostellar embryo forms following the 
non-homologous collapse of a much larger gravitationally 
unstable core, and therefore it accretes mainly from its own 
co-moving placenta. At the other extreme, it is argued that 
a protostellar embryo quickly becomes sufficiently decou- 
pled from the ambient gas that it can roam around compet- 
ing with other embryos for the same reservoir of accretable 
gas. For the purpose of this section, we assume that nature 
cleaves to the second extreme. 

7.1. Competitive accretion in gas-rich proto-clusters 

Once an ensemble of protostellar embryos has formed, 
still deeply embedded in its parental prestellar core and/or 
parental disc, the individual embryos evolve by accreting 
gas, and by interacting dynamically with one another. The 
accretion histories of individual embryos differ due to their 
varying circumstances, leading to a spectrum of protostel- 
lar masses, extending from high masses down to below the 
H-burning limit. Those which spend a long time moving 
slowly through the dense gas near the centre of the core 
can grow to high mass. Conversely, those that spend most 
of their time moving rapidly through the diffuse gas in the 
outer reaches of the core do not grow much. This process 
of 'competitive accretion' may be a major element in the 
origin of the IMF, as first pointed out by Zinnecker (1982). 

Over the past decade, many numerical simulations of the 
formation of star clusters by fragmentation and competi- 
tive accretion have been performed (e.g. Chapman et al, 
1992; Turner et al, 1995; Whitworth et al., 1995; Bonnell 
et al., 1997, 2001a, 2001b; Klessen et al, 1998; Klessen 
andBurkert, 2000, 2001; Bate et al, 2002a, 2002b, 2003; 
Goodwin et al. 2004b, 2004c). Bonnell et al. (1997, 2001a, 
2001b) have shown through numerical simulations and an- 
alytical arguments that competitive accretion in large-iV 
proto-clusters can reproduce the general form of the IMF. 

7.2. Unstable multiple systems 

Multiplicity studies of main sequence and evolved stars 
have revealed that 15 to 25% of all stars, when studied in 
sufficient detail, are triple or higher-order multiples (e.g. 
Tokovinin and Smekhov, 2002; Tokovinin, 2004). It follows 
that the formation of multiple stars is an important element 
of the star formation process. 

The multiplicity fraction among PMS stars is poorly 
known, partly because of the difficulty in studying stars in 



the embedded phase, but it appears to be at least as high as 
- and probably much higher than - for more evolved stars 
(e.g. Reipurth, 2000; Looney et al, 2000; Koresko, 2002; 
Reipurth et al, 2002, 2004; Haisch et al, 2004), see also 
the chapters by Duchene et al, Goodwin et al. and Bur- 
passer et al. in this volume. Most stars are formed in em- 
bedded clusters, and there is increasing evidence that the 
primary building blocks of clusters are small subclusters of 
< 20 stars which quickly dissolve and merge to form the 
more extended cluster (e.g. Teixeira etal, 2006). 

Nonhierarchical multiple systems, in which the time- 
averaged distances between components are comparable, 
are inherently unstable (e.g van Albada, 1968). Within 
about a hundred crossing times a triple system is likely 
to have ejected one member, most likely the least massive 
component, since the ejection probability scales approxi- 
mately as the inverse third power of the mass (e.g. Anosova, 
1986; Mikkola and Valtonen, 1986). Although most non- 
hierarchical systems disintegrate in this way, the existence 
of numerous stable hierarchical triple systems shows that 
this is not always the case. Ejected members leave with 
a velocity that, to first order, is comparable to the veloc- 
ity attained at pericenter in the close triple encounter, and 
depends on the geometry of the encounter, the energy and 
angular momentum of the system, and the masses of the 
components (e.g. Standish, 1972; Monaghan, 1976; Sterzik 
and Durisen, 1995, 1998; Armitage and Clarke, 1997). 

The disintegration of a small multiple system is most 
likely to occur during the deeply embedded Class phase, 
while massive accretion from a surrounding envelope is still 
taking place (Reipurth 2000). If the ejection leads to an es- 
cape, then the accretion halts and the final mass of the object 
is capped {Klessen, Burkert, & Bate 1998). 

7.3. Dynamical ejection as a source of brown dwarfs 

If a protostellar embryo is ejected from its natal core with 
a mass below the H-burning limit, then it becomes a brown 
dwarf (Reipurth and Clarke, 2001). All that is needed for 
this to happen is for a prestellar core to spawn more than 
two protostellar embryos; for at least one of them to be less 
massive than 0.075 M at the outset; and for one of them 
to stay less than 0.075 M long enough to be ejected by 
dynamical interaction with the other embryos. 

Simulations suggest that forming more than two proto- 
stellar embryos in a collapsing core is routine, as is the ejec- 
tion of brown dwarfs and very low-mass stars (e.g. Bate 
et al, 2002a; Delgado-Donate et al, 2003, 2004; Good- 
win et al, 2004a, 2004b). However, these simulations may 
be misleading, (i) They use sink particles (Bate, Bonnell 
and Price, 1995), and thereby create protostellar embryos 
which at the outset are inevitably low-mass (~ 0.005 M ), 
very prone to dynamical interaction (being effectively point 
masses, albeit with gravity softening), and unable to merge. 
Thus, although formation by ejection seems inevitable, it 
may be less efficient than these simulations suggest. As 
discussed by Goodwin and Kroupa (2005) and Hubber and 
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Whitworth (2005), the observed binary statistics are incom- 
patible with too many ejections, and the number of proto- 
stellar embryos undergoing dynamical interactions within a 
single prestellar core should be relatively small {J\f < 4). 
(ii) They do not include magnetic fields. Boss (2002) ar- 
gues that a magnetic field may promote fragmentation of 
a collapsing core, by inhibiting the formation of a central 
density peak; however, his code does not capture all the 
possible MHD effects, in particular the anisotropy of mag- 
netic pressure and the torques exerted by a twisted field. In 
contrast Hosking and Whitworth (2004) have simulated the 
collapse of a rotating core with imperfect MHD, and find 
that fragmentation is inhibited. More work is needed. 
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Fig. 2. — Four stars formed in a low-mass core (~ 5 M Q ) with a 
low initial level of turbulence (-E TURB ~ 0.05 |-B GRAV |) {Good- 
win et al. 2004b). Near the centre is a triple system containing 
a close binary (barely resolved pair of black dots). The object 
towards the upper righthand corner is a brown dwarf which has 
been ejected and is unbound from the triple. It has a significant 
disc (M DISC ~ 0.001 M ). The frame is 10 3 AU across, and the 
time is 0.073 Myr since the start of collapse. 

In their original paper, Reipurth and Clarke (2001) con- 
jectured that brown dwarfs formed by ejection could have 
a higher velocity dispersion than more massive H-burning 
stars, and that this might be detected, either in observa- 
tions of the velocity dispersions of young star clusters, or 
as a diaspora of brown dwarfs around more evolved clus- 
ters. They also pointed out that violent ejections would re- 
sult in smaller accretion discs, and therefore a shortened 
quasi-T Tauri phase. Numerical simulations of star clus- 
ter formation in highly turbulent massive cores (e.g. Bate 
et ah, 2003; Bate and Bonnell, 2005) find that the brown 
dwarfs do indeed have quite a high velocity dispersion 
(~ 2 to 4kms~ 1 ), but they are difficult to distinguish from 
the low-mass H-burning stars, because both are frequently 



ejected from small dynamically-unstable groups. Simula- 
tions of low-density star-forming regions by Goodwin et 
al. (2005) find that the velocity dispersions are typically 
somewhat lower (~ 0.5 to lkms -1 ), and again the ve- 
locity dispersions of brown dwarfs and H-burning stars are 
hard to distinguish, because they are partially masked by 
the velocity dispersion between the different cores in which 
small subgroups (< 5) of stars are born. For the same 
reasons, segregation of brown dwarfs from H-burning stars 
may be difficult to detect. Moreover, given their low ejec- 
tion speeds, some brown dwarfs formed by ejection retain 
significant discs, and are therefore presumably well able 
to sustain accretion and outflows, as observed (see §2.3 
and the chapter by Luhman et al.). Thus kinematical and 
spatial information on brown dwarfs, and the signatures 
of accretion and outflow, provide less powerful constraints 
on brown dwarf formation than were initially surmised by 
Reipurth and Clarke (2001). 

7.4. Dynamical interactions in early stellar evolution 

The binary statistics of brown dwarfs provide another 
potential constraint on formation mechanisms (see §2.2 and 
the chapters by Burgasser et al., Duchene et al., Goodwin 
et al, and Luhman et al.). Formation by ejection may be in- 
compatible with the relatively high frequency of close BD- 
BD binaries, if this is confirmed (e.g. Maxted and Jeffries, 
2005; Joergens, 2005a). As well as resulting in the ejec- 
tion of some brown dwarfs, dynamical interactions must 
also influence the binary statistics (distributions of primary 
mass, mass ratio, semi-major axis and eccentricity) of the 
brown dwarfs and low-mass H-burning stars that do not get 
ejected (e.g. Sterzik and Durisen, 2003; Kroupa and Bou- 
vier, 2003). Hubber and Whitworth (2005) have shown that 
if they take the observed distributions of core mass, radius 
and rotation rate, convert each core into a ring of 4 to 5 
stars with masses drawn from a log-normal distribution hav- 
ing dispersion cr logio[M] = 0.6, and then follow the dissolu- 
tion of the ring by pure A^-body dynamics, they reproduce 
rather well the observed distribution of multiplicity and bi- 
nary statistics in young clusters, as a function of primary 
mass. Umbreit et al. (2005) have investigated the disinte- 
gration of nonhierarchical accreting triple systems, and find 
that they are able to produce the observed separation dis- 
tribution of close binary brown dwarfs. Thus dynamical 
interactions may make important contributions, both to the 
formation of brown dwarfs, and to their binary statistics. 

In addition to producing brown dwarfs and very low 
mass stars, and helping to shape the lower-mass end of the 
IMF and the statistics of binary and higher multiple sys- 
tems, dynamical interactions may be the key to understand- 
ing a variety of other phenomena in early stellar evolution. 
Reipurth (2000) notes that the different sizes and separa- 
tions of shocks in Herbig-Haro flows can be understood as a 
fossil record of the dynamical evolution of a newly formed 
binary. He also notes that dynamical interactions may on 
occasion lead to a departure from the standard evolutionary 
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picture of a star passing smoothly from the Class stage 
through the Class III stage. Instead stochastic dynamical in- 
teractions may lead to the sudden ejection of an object from 
the Class or I stage, resulting in its abrupt appearance as 
a Class II or even a Class III object. The infrequent young 
binaries with infrared companions may be related to such 
events. Finally, the FU Orionis eruptions may be related 
to the formation of a close binary, and could result from 
the viscous interactions of circumstellar material as the two 
components in a newly formed binary spiral together (Bon- 
nell and Bastien, 1992; Reipurth and Asp in, 2004). 

8. PHOTO-EROSION OF PRE-EXISTING CORES 

A fifth - and somewhat separate - mechanism for form- 
ing brown dwarfs is to start with a pre-existing core of stan- 
dard mass (i.e. < M ) and have it overrun by an HII re- 
gion (Hester et ai, 1996). As a result, an ionisation front 
(IF) starts to eat into the core, 'photo-eroding' it. The IF 
is preceded by a compression wave (CW), and when the 
CW reaches the centre, a protostar is created, which then 
grows by accretion. At the same time, an expansion wave 
(EW) is reflected and propagates outwards, setting up the 
inflow which feeds accretion onto the central protostar. The 
outward propagating EW soon meets the inward propagat- 
ing IF, and shortly thereafter the IF finds itself ionising gas 
which is so tightly bound to the protostar that it cannot be 
unbound by the act of ionisation. All the material interior 
to the IF at this juncture ends up in the protostar. On the 
basis of a simple semi-analytic treatment, Whitworth and 
Zinnecker (2004) show that the final mass is given by 

~0 01M ( * ) 6 ( Y ~ 3 ( n a y h 
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(30) 

where a, is the sound speed in the neutral gas of the core, 
A/" LyC is the rate at which the exciting star(s) emit ionising 
photons, and n Q is the density in the HII region. 

This mechanism is rather robust, in the sense that it pro- 
duces very low-mass stars for a wide range of initial con- 
ditions, and these conditions are likely to be realised in 
nature. Indeed, the evaporating gaseous globules (EGGs) 
identified in M16 by Hester et al. (1996) - and subse- 
quently in other HII regions - would appear to be pre- 
existing cores being photo-eroded in the manner we have 
described (e.g. McCaughrean & Anderson, 2002). How- 
ever, the mechanism is also very inefficient, in the sense 
that it usually takes a rather massive pre-existing prestel- 
lar core to form a single brown dwarf or very low-mass H- 
burning star. Moreover, the mechanism can only work in 
the immediate vicinity of an OB star, so it cannot explain 
the formation of all brown dwarfs, and another mechanism 
is required to explain those seen in star formation regions 
like Taurus. Nonetheless, if the majority of stars are born 
in Trapezium like clusters, rather than Taurus-like regions, 
then photo-erosion should remain a contender for produc- 
ing some brown dwarfs. Brown dwarfs formed by photo- 



erosion should include close BD-BD binaries. It is unclear 
whether they can retain significant accretion discs. 

9. SIMULATIONS OF CLUSTER FORMATION 

With several different likely mechanisms for the produc- 
tion of brown dwarfs, the question arises: which, if any, is 
the dominant formation mechanism? This is only likely to 
be answered through numerical simulations that are able to 
model the full star formation process, including all the rel- 
evant physical ingredients (gravity, hydrodynamics, mag- 
netic fields, radiative transfer, and chemistry). There is a 
huge effort underway to perform such simulations, but there 
are formidable numerical challenges to overcome. 

9.1. Turbulent cloud collapse 

The most comprehensive simulations to date are those 
of Bate, Bonnell and Bromm (2002a, 2003), Bate and Bon- 
nell (2005) and Bate (2005). These model the collapse 
and fragmentation of turbulent molecular clouds with mass 
M ~ 50 M , initial diameter 0.18 to 0.38 pc and initial 
temperature T ~ 10 K, to form small stellar clusters con- 
taining ~ 50 stars, including numerous brown dwarfs. The 
clouds are seeded with a power spectrum of supersonic ve- 
locity structure that matches the scaling of velocity disper- 
sion with length-scale observed in molecular clouds (Lar- 
son, 1981) and is allowed to decay during the simulations. 
The key difference between these simulations and earlier 
ones is that they are able to produce large numbers of ob- 
jects from which statistical quantities can be derived (e.g. 
the form of the IMF), but simultaneously they also resolve 
down to the opacity limit for fragmentation (Section 3) and 
so they are able to follow the formation of all the stars 
that the clouds produce. They also resolve gaseous discs 
with radii down to ~ 10 AU and binaries with separations 
greater than ~ 1 AU. On the other hand, they do not include 
magnetic fields, radiative transfer, chemistry or feedback. 
Therefore, for instance, they are unable to investigate the 
fraction of brown dwarfs that might form via photo-erosion. 

With these caveats in mind, the simulations generate 
clusters which are quite realistic. Starting with initial con- 
ditions typical of the molecular clouds in our Galaxy, the 
simulations produce roughly equal numbers of stars and 
brown dwarfs, with an IMF which is roughly compatible 
with that observed, at low masses; the high-mass end of 
the IMF (> M ) is not usefully constrained by these sim- 
ulations. All stars, including brown dwarfs, originate from 
the fragmentation of dense filaments of molecular gas and 
from disc fragmentation, as illustrated in Fig. 3. Crucially, 
however, those that end up as brown dwarfs are those that 
avoid accreting large amounts of gas from the surround- 
ing cloud. They are able to avoid becoming stars because 
they are ejected dynamically from unstable multiple sys- 
tems, thereby terminating their accretion. In general it is 
easier to form brown dwarfs by disc fragmentation, because 
the resulting protostellar embryos are then born in a dense 
multiple system, and they can therefore be ejected quickly, 



14 





3 




B /^2 * 


A 


A 


A 






A 


.* 







Fig. 3. — Two sequences illustrating the brown dwarf formation mechanisms that occur in the simulations of Bate et al. (2002a, 2002b, 
2003). The upper sequence shows two brown dwarfs (square and triangle) forming in a circumbinary disc, while the lower sequence 
shows two brown dwarfs (square and triangle) forming in a filament. In both cases, these objects remain as brown dwarfs because they 
are dynamically ejected before they can accrete sufficient mass to ignite H-burning. In the upper case, they are ejected from a multiple 
system formed by disc fragmentation. In the lower case, the two brown dwarfs form in separate filaments, but then fall into, and are 
ejected from, the multiple system that exists at the intersection of the two filaments. Each panel is 600 AU across. 



before acquiring too much mass. In contrast, protostellar 
embryos formed by filament fragmentation are born in rel- 
ative isolation, and they first have to fall down the filament 
into a dense cluster before they are ejected; their final mass 
is therefore less likely to stay below the H-burning limit. 
Thus, these simulations support the ejection hypothesis for 
the origin of brown dwarfs, with disc fragmentation also 
playing an important role in generating dynamically unsta- 
ble systems. 

The brown dwarfs formed in these simulations do not 
frequently have companions. With three individual calcula- 
tions now published, the overall frequency of BD-BD and 
very low-mass binaries is ~ 5%. Most of these systems are 
closer than ~ 20 AU. Sun-like stars with brown dwarf com- 
panions at separations less than 10 AU are also very rare, 
at ~ 2%. Wider brown dwarf companions at hundreds or 
thousands of AU are much more common, although none of 
these systems have reached dynamic stability when the sim- 
ulations are terminated. Most of these results are consistent 
with current observations, with the exception of the BD- 
BD binary frequency, where the value from the simulations 
is about three times lower than the observed value (~ 15%). 
There are at least two possible reasons for this. First, most 
observed and calculated BD-BD binaries have separations 

< 20 AU, but the simulations are unable to resolve discs 

< 10 AU, so the discrepancy might be solved with better 
resolution. If this is not the case, missing physics may 
be the problem (e.g. the effects of radiative transfer on 
disc fragmentation). Finally, only one brown dwarf that 
was ejected during any of the simulations had a resolved 
disc (radius R > 10 AU). This implies that brown dwarfs 
formed from highly turbulent cores should only have small 



discs. In contrast, the simulations of cores with low levels 
of turbulence performed by Goodwin et al. (2004a, 2004b) 
produce brown dwarfs with somewhat larger discs (see Fig. 
2). This indicates that disc size may be a function of birth 
environment. Although many brown dwarfs are observed 
to have discs from their spectral energy distributions, the 
distribution of their sizes is not currently known. 

10. CONCLUSIONS 

Since the statistical properties of brown dwarfs appear to 
form a continuum with those of low-mass H-burning stars 
(§2), we have argued that brown dwarfs form like stars, that 
is to say, on a dynamical timescale and by gravitational in- 
stability, with an homogeneous initial elemental composi- 
tion, the same as the interstellar medium from which they 
form. In this regard brown dwarfs are distinct from plan- 
ets, which we define to be objects that form on a longer 
timescale, by the accumulation of a rocky core and - if cir- 
cumstances allow - the subsequent acquisition of a gaseous 
envelope, leading to an initially fractionated elemental com- 
position and a deficit of light elements. 

We have evaluated the minimum mass for a brown 
dwarf (§3) by considering the cooling required for a very 
low-mass prestellar core to condense out on a dynamical 
timescale. We have treated several different scenarios: hier- 
archical 3-D fragmentation; one-shot 2-D fragmentation of 
a shock-compressed layer; and fragmentation of a Toomre- 
unstable disc. All three cases yield values of M MIN in the 
range 0.001 to 0.004 M for contemporary star formation 
in the solar vicinity. This suggests that there may be some 
overlap between the range of masses occupied by stars and 
planets. In hotter environments and at earlier epochs, M MIN 
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was probably larger, and brown dwarfs were therefore less 
common. We also suggest that the thermodynamics of discs 
make it easier for proto brown dwarfs to condense out at 
large radii > 100 AU, and that these proto brown dwarfs 
may fragment to produce close BD-BD binaries, due to the 
dissociation of H 2 and the opacity gap at T ~ 2000 K. 

We have discussed five possible mechanisms for form- 
ing brown dwarfs. Turbulent fragmentation of molecular 
clouds may deliver prestellar cores of such low mass that 
their subsequent collapse (and possible fragmentation) can 
only yield brown dwarfs (§4). The collapse and fragmen- 
tation of more massive cores is likely to deliver protostel- 
lar embryos with a wide range of masses, many of which 
are too low to support hydrogen burning (§5). Similarly, 
disc fragmentation is likely to deliver low-mass protostel- 
lar embryos (§6). These protostellar embryos may undergo 
competitive accretion (as a result of which some of them 
evolve to higher mass) and dynamical interactions (as a re- 
sult of which some of them are ejected before they reach 
H-burning masses) (§7). Lastly, cores which find them- 
selves overrun by an HII region may be photo-eroded by 
the resulting ionisation front and end up spawning brown 
dwarfs (§8). None of these mechanisms is mutually exclu- 
sive, and in the most advanced simulations of cluster for- 
mation (§9) collapse and fragmentation, disc fragmentation, 
competitive accretion and dynamical ejection all occur con- 
currently. 

However, even these simulations do not capture all the 
deterministic physics. It requires a fully radiative, effec- 
tively inviscid, 3-D magneto-hydrodynamical simulation to 
evaluate properly the thermal effects which influence the 
minimum mass for star formation, the angular momentum 
transport processes which influence the binary statistics, 
and the TV-body dynamics which influence the clustering 
properties of brown dwarfs. Work to develop and validate 
such codes is ongoing. 

We believe that all the proposed mechanisms operate in 
nature, and that once they have been properly modelled, the 
ultimate task will be to determine their relative contribu- 
tions to the overall brown dwarf population. These relative 
contributions may depend on environment, metallicity and 
epoch, and may therefore lead to local and/or cosmological 
variations in the ratio of brown dwarfs to H-burning stars, 
and in the binary and accretion statistics of brown dwarfs. 
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